Goos-Hanchen and Imbert-Fedorov Shifts of Vortex Beams at Air-Left-Handed 

Material Interfaces 
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In this paper, we present a systematic study of beam shifts and angular momenta of paraxial 
vortex beams at air-left-handed material (LHM) interfaces. It is shown that, compared to their 
counterparts at air-right-handed material (RHM) interfaces, the spatial Goos-Hanchen (GH) and 
Imbert-Fedorov (IF) shifts remains the same, while the angular GH and IF shifts are reversed, at air- 
LHM interfaces. The spatial and angular shifts of paraxial vortex beams have their respective origins 
in transverse angular momenta and transverse linear momenta. The spatial GH and IF shifts remain 
unreversed as a result of the both reversions of transverse angular momenta and z-component linear 
momentum, while the angular GH and IF shifts are reversed due to that the z-component linear 
momentum is reversed and the transverse linear momenta are unreversed, at air-LHM interfaces. In 
addition, we perform a quantitative analysis on spin-orbit angular momentum conversion and orbit- 
orbit angular momentum conversion, which further helps us understand the essence of vortex beam 
shifts at air-LHM interfaces and their fundamental distinctions with those at air-RHM interfaces. 

PACS numbers: 42.25.-p, 41.20.Jb, 78.20.Ci 



I. INTRODUCTION 

The reflection and transmission of light at air-medium 
interface has been a topic of ongoing investigation and 
contention. Initially, the amplitudes and directions of re- 
flected and transmitted light were determined by Fresnel 
formulas and Snell's law. Deviations were later found 
in total internal reflection (T'lR), however. They are 
Goos-Hanchen (GH) shift [lHi| and Imbert-Fedorov (IF) 
shift (5|-[l0|, which are parallel and perpendicular to the 
incident plane, respectively. Artmann formula Q for GH 
shift and Schilling formula @ for IF shift in TIR are 
widely accepted. The exact formulas for IF shifts in par- 
tial reflection and transmission, however, have divided 
many physicists [Tl| - [T^ . Owing to Hosten's precise 
measurement, the debate was settled down jl5l] . Apart 
from the constant spatial GH and IF shifts, the light 
beam also experiences angular shifts [l6l,fl7|, which in- 
crease proportionally with propagation distance z. In 
general, the spatial GH and IF shifts stem from spin- 
orbit interaction at the air- medium interface [Tsl . [l9j . To 
satisfy the z-componcnt angular momentum conservation 
law, the reflected and transmitted beams must possess 
extrinsic transverse orbital angular momenta, which re- 
sults in spatial GH and IF shifts. The angular shift is es- 
sentially a difFractive correction on light beams and gov- 
erned by the linear momentum conservation law [T8l. l20||. 

The emergence of left-handed materials (LHM) has 
brought about great opportunities and sophisticated 
pathways to manipulate light [21-24j. As a branch of 
metamaterials, LHM usually refers to those materials 
whose permittivity and permeability are negative simul- 
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taneously. It has shown very unique properties, like 
negative refraction [2^, inverse Doppler effect [1^ [26| . 
unreversed rotational Doppler effect [13] and inverse 
Cherenkov radiation [28l ). Apart from these properties, 
beam shifts in LHM have been demonstrated to be signif- 
icantly different from right-handed materials (RHM) |29l - 
|33| . Theoretical papers [1^ [sOj indicate that light beam 
experiences negative GH shift at air-LHM interface and 
this shift can be amplified in layered structure. The 
IF shift in partial refiection, however, remains unre- 
versed owing to unreversed spin angular momentum in 
LHM [33*] . As for vortex beam, it carries intrinsic orbital 
angular momentum. Therefore, orbit-orbit conversion is 
inevitable in refiection and transmission. Will the spa- 
tial GH and IF shifts remain unreversed in the presence 
of orbital angular momentum? How about the angular 
shifts? What does the physical picture of spin-orbit and 
orbit-orbit conversions look like? Clarifying these prob- 
lems is not only conducive to understanding the detailed 
behaviors of vortex beams at air-LHM interface but also 
sheds a new light on the essence of angular momentum 
and linear momentum of vortex beam. 

In this paper, we endeavor to offer concrete expressions 
of spatial and angular shifts of vortex beam at air-LHM 
interface and elaborate on their relations with the lin- 
ear and angular momenta. We also contrast these results 
with air-RHM interface. Although some of the results are 
similar to our previous paper [33|, we conduct consider- 
ably rigorous analyses on the topics and thorough expla- 
nations on the results. The rest of the paper is arranged 
as follows. In SecHH we adopt angular spectrum method 
to derive the electric fields of reflected and transmitted 
vortex beams at air-LHM interface. Except for some spe- 
cial cases, for instance, in the vicinity of Critical angle of 
TIR and Brewster angle, the expressions of electric fields 
generally hold true. The longitudinal fields are included 
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as well, which take on polarization-sensitive vortex struc- 
ture l56| , providing a new perspective on the mechanism 
of IF shifts other than spin-orbit conversion. In Sec. IIIH 
we adopt the operator method to calculate the spatial 
and angular shifts, which is significantly different from 
conventional calculation method 113, HB] and save 

us considerable time and efforts. The GH and IF shifts 
in TIR and partial transmission are demonstrated and 
contrasted with air-RHM interface. The impact of inci- 
dent angle on beam shifts is analyzed as well. In Sec. lIVi 
we calculate the linear and angular momenta of incident, 
reflected, and transmitted beams, respectively. The ori- 
gins of spatial and angular shifts are clarified. The quan- 
titative analysis on spin-orbit conversion and orbit-orbit 
conversion is also demonstrated. 




II. ELECTRIC FIELDS OF REFLECTED AND 
TRANSMITTED VORTEX BEAMS 



In this study, we adopt the angular spectrum method 
to establish a model for the reflected and transmitted 
vortex beams at air-LHM interface. The main proce- 
dures of this method are as follows. We first decompose 
the vortex beams into plane waves with finite spectral 
width. Then, we analyze the incident angle and Fres- 
nel coefficients of each plane wave component separately. 
In this case, the incident angle and Fresnel coefficients 
of each plane wave slightly differ from that of the main 
Fourier component (also known as central wave compo- 
nent). Therefore, we expand the Fresnel coefficients in 
Taylor series around the central incident angle and make 
approximation to the first order. Afterwards, we trans- 
form the electric fields from momentum space to position 
space. 

The geometry of reflection and transmission is demon- 
strated in Fig. [TJ The incident vortex beam propagates 
along Zi axis and impinges on the air-LHM interface. The 
incident plane of the main Fourier component is xoz. The 
reflected wave and transmitted wave of the main Fourier 
component travel along the axis and zt axis, respec- 
tively. The incident angle, reflection angle, and trans- 
mission angle of the main Fourier component are Oi, 6r, 
and 9t, respectively. The relations among the coordi- 
nates O - Xyz, Oi - XiyiZi, Or - XrUrZr, and Ot - XtVtZt 

are determined by the Snell's Law: 
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where r = i, r, or [x y z] and [x^. y^- z^ ] are 
the unit basis vectors of the coordinates o — xyz and Or — 
Xrj/rZr, respectively, dr '&t]^[di tt - 6i -Ot], 
9i = Or and sin^t — sin0i/|n|. Constant n is the re- 
fractive index of the LHM. For an arbitrary plane wave, 
we assume that the incident plane is XOZ. The inci- 
dent, reflected, and transmitted plane wave travel along 



FIG. 1. (Color online) Geometry of reflection and transmis- 
sion. The incident plane of the main Fourier component is 
xoz (marked with green lines). For an arbitrary plane wave, 
the incident plane is XOZ (marked with blue lines). 



the Zi, Zr, and Zt axis, respectively. The incident an- 
gle, reflection angle, and transmission an^le of an arbi- 
trary Fourier component are 0^, 0^, and 0^, respectively. 
The relations among coordinates O — XYZ, Oi — XiYiZi, 
Or — XrYrZr, and Ot — XtYtZt are easily acquired by 
revising [^i §r M in Eq. ©as [o'. tt - 0[ -O't]. 
From Fig. [1] we derive the incident angle of an arbitrary 
plane wave, Oi = arccos (— fc^;. sin^i/A: -I- fc^. cos^i/fc), 
where k is the wave number in vacuum, k^^, ky., kz^ are 
the components of wave vector along the Xi, yi, Zi axis, 
respectively. We can expand the incident angle Oi around 
the central incident angle Oi in series of kx^/k and kyjk. 
Therefore, 9^ ~ Oi + k^^jk, Ot ^ Ot + kxj{\n\r]k), where 
7] = cos Ot I cos Oi . 

After introducing the geometry of reflection and trans- 
mission, we start analyzing the angular spectrum of re- 
flected and transmitted beams by using transformation 
matrix. Note that the incident, reflected, and transmit- 
ted vortex beams are presented in local coordinate sys- 
tems 0-j- X -J- y -J- z -J- . The angular spectrum of incident beam 
is: 



axi + /3y, - 1 {akx, + (iky^ ) z^ 



2 

X exp [ 



wq {-ikx, + sgii[l]ky.) / V2 
{kl+kl)/4]. 



Ui, 



(2) 



where Ui is the angular spectrum of vortex beams, / is 
the vortex charge, C; = ^/2J{^^\l\l) is the normalization 
constant, wq is the width of beam waist, k^^ and ky. 
are wave vector components along the Xi and yi axis, k 
is wave number in vacuum, sgn[/] is the sign function, 
a and (3 are Jones vectors, |ap -I- |/3p = 1. There are 
two parameters that characterize the polarization state of 
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paraxial beams: cr = 2 Im [a* /3] , x = 2 Re [a*/3]. a is the 
degree of circular polarization. A value of cr = +1 cor- 
responds to left circularly polarized light beam, whereas 
a parameter of cr = —1 stands for right circularly polar- 
ized light beam, cr = represents linear polarization and 
values between and 1 should correspond to elliptically 
polarized states, x is the degree of linear polarization. It 
is generally recognized that the circularly polarized vor- 
tex beam has intrinsic angular momentum {I + <j)Ti per 
photon [H,!!!. 

We first write the electric field in coordinate Oi — 
XiYiZi. This target would be achieved in three steps. 
First, we write the electric field in o — xyz by using the 
matrix in Eq. ([T]). Second, we transform the electric field 
from o — xyz to O — XYZ . This step is accomplished by 
using the relations: X = (zxki) x z/| (zxki) xz|, Y = zx 
ki/|z X ki|, Z = z, where kj = (fc^iX^ -I- ky^yi + k^^z^) /k. 
Therefore, this transformation matrix is given below: 



x 




Y 




Z 





i~\-kz ■ sin 6 



k sin 6 

k-u. 



k sin 9 ■ 

cos Oi-{~kz ■ sin 6 







Third, 
to 0^ 



k sin 6 ■ k sin 9 . 

' ' 

(3) 

we transform the electric field from O — XY Z 
- XiYiZi using revised version of Eq. ([T]). The 



transformation matrix of these three steps is written as 



1 

■ cot i 



cot Bi 
~k 
1 



(4) 



In Eq. (HI), we assume that terms equivalent to or higher 
than {kxjk)"^, {ky./k)"^ are null. 

After we get the expression of electric field in Oi — 
XiYiZi coordinate, the refiected and transmitted fields 
can be easily obtained by multiplying the Fresnel coef- 
ficients. We define the amplitude reflection coefficients 
of the main Fourier component as r^, r^, the ampli- 
tude transmission coefficients of the main Fourier com- 
ponent as tp, ts, where p and s denote p-polarized and 
s-polarized state. We also define the amplitude reflec- 
tion coefficients of an arbitrary plane wave as Tp, r^, the 
amplitude transmission coefficients of an arbitrary plane 
wave as t^, t^. To simplify the Fresnel coefficients of an 
arbitrary wave component, we expand them around the 
central incident angle 6i in series of kx^/k and ky./k, and 
retain the first order term. Thus, r„ 



Note that, the reflected and trans- 



dOi k 

mitted electric fields are still presented in Or 
and Ot — XtYtZt coordinates. Hence, we need to trans- 
form them back to the Or — XrVrZr and Ot — XtytZt coor- 
dinates. The transformation matrixes are 
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kjj . cot (. 



ky . cot (. 
^~k 
1 



yt 



1 

. cot Si 



k 

k^i 

717] k 



-rjky . cot 9i 
% 

1 

nk 



nr}k 



1 





X, 




Y, 




z. 



(5) 



As 



The detailed procedures are similar to the Eq. ^ 
mentioned previously, we eliminate terms equal to or 
higher than (kxi/k)'^ and (kyjk)'^. 

The Fresnel coefficients only reveal the amplitudes of 
reflected and transmitted vortex beams. To get the full 
expressions of reflected and transmitted vortex beams, 
we still need the phase matching conditions. For an 
arbitrary plane wave component, the boundary condi- 
tions [36] require that in the plane xoy, expfzkj • r] = 
exp[ik,. -r] = exp[ikj-r], where k;, k^, and kt are the wave 
vectors of incident, reflected, and transmitted beams, re- 
spectively. Therefore, the phase matching conditions are 
presented like this: 

— — Ixi — ■ 

kxt ~ ^Xi/Vi ^Vt ~ ^yn ~ '^^Zi ■ (6) 

Combining the Fresnel coefficients and Eqs. (0) and 
([6]) , we get the electric fields of reflected and transmitted 
beams: 



a r, 



(3 r. 



Wi^ 



+ /3 {rs + rp) cot 0i-^ 



a [rs + Tp) cot 9i 



- {aVpkx^ + Prsky^) UrZr, 



a(ip + ry||%)+/3fe-ryt.)cot0.^ 
^(i. +?7||^j +a{7jtp-ts)cote,'^ 



ik 



{atpkxt + ptskyJutZt, 



Urfr 

(7) 

utYt 
(8) 



where Ur,t = {lr,tkx,.^t, ky^,t) ^ 7r = -1, 7t = V- 
Note that the reflected and transmitted vortex 
beams experience extra phase shifts [13] associ- 
ated with diffraction while propagating, which are 



exp 



1 - 



and exp 



inkzt 1 



respectively. 

At this stage, we feel obliged to deliver several com- 
ments. First, the phase matching conditions, as demon- 
strated in Eq. ([6]), play a crucial role in this paper. For 
instance, kx^ = —kx^ and ky^ = ky. show that the 
wavefront of reflected vortex beam is reversed (shown 
in Fig. [2|), resulting in a reversed orbital angular mo- 
mentum —Ih per photon, kz^ = nk^. shows that the 
transmitted beam in LHM undergoes negative phase ve- 
locity, resulting in a reversed helical wavefront (shown 
in Fig. [5]). The reversed wavefront in LHM does not 
lead to reversed orbital angular momentum, which is ex- 
plained in Sec. IIVI The equation kx^ = kxi/rj means the 
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FIG. 2. (Color online) The helical wavefronts of incident, 
reflected, and transmitted beams. Incident, reflected, and 
transmitted vortex beams propagate along Zi, Zr, Zt axis, re- 
spectively. We choose I — 1 for incident vortex beam. The 
incident, reflected, and transmitted vortex beams have anti- 
clockwise, clockwise, and clockwise helical wavefronts, respec- 
tively. 



transmitted beam experiences a beam deformation in the 
cross section. Second, in the vicinity of Critical angle 
Oc of TIR and Brewster angle 63 {\0, - 6Ib,c| < 6*0, 
where Bq — i/ikw^)) is the beam divergence angle), we 
should be cautious. For example, in the vicinity of Crit- 



ical angle \9i 



< 60), the first derivatives of 



are extremely large 38], meaning that Eqs. ([7]) and 
should be revised to include terms higher than {kxJkY 
and {ky./k)'^. If a p-polarized vortex beam impinges on 
the air-LHM interface with incident angle being in the 
vicinity of Brewster angle {\0i — 9b\ < 6*0), the power of 
reflected beam would be in the scale of 9q (approximately 
10~* ~ 10~^ for paraxial beams) comparing with the in- 
cident beam. This constitutes the reason why the exper- 
iments [39', '4^ fail to collect reliable data concerning the 
shift of p-polarized beams in the vicinity of Brewster an- 
gle. Third, although the beam shape and shift of reflected 
vortex beam are already clear [3; Esf . The profile of 
transmitted vortex beam is far from clear-cut. We formu- 
late a rigorous transmitted field in the appendix of our 
paper. By referencing this appendix, the beam shape of 
the transmitted vortex beam is easily acquired. 



III. SPATIAL AND ANGULAR SHIFTS 

The paraxial wave equation is identical with two- 
dimensional Schrodinger equation with z replaced b y t . 
Hence, the operator formalism proposed by Stoler 1411 
has become a powerful tool in physical optics [42H44I ] . 



In this section, we will apply the operator formalism to 
calculate the centroid of the reflected and transmitted 
beams. In momentum space, the transverse position op- 
erator is idk^ . Therefore, the centroid of beams in a 
given plane z = const is readily given by 



(GE|GE) ' 



(9) 



where r_L = xx + yy, = gl^x -I- g|-y, the prop- 
agation operator G — exp [inkz — (k^. -\- /c^)] . The 
above equation can be easily formulated into the follow- 
ing one: 



(EK^k^lE) 
(EIE) 



z (E|k^|E) 
nk (EIE) 



(10) 



where k is the wave number in vacuum and n is the re- 
fractive index. The first and second terms of Eq. pUj) are 
the spatial and angular shifts, which are independent and 
dependent on z, respectively. From the second term of 
Eq. (fTO|) . we can easily find that the negative refractive 
index n results in reversed angular shifts in LHM [33{ . 
For the reflected and transmitted beams, the GH and IF 
shifts can be expressed as: {xt) ~ Aa^r + Zt^Ox^ and 
{yr) = Ayr + ZrA6y^ , where t = r,t. 

We first calculate the spatial and angular shifts of re- 
fiected beams. Substituting Eq. ([T]) into Eq. (|10l) . we get 
the spatial and angular shifts of reflected beams 



AXr 



xl cot 9i (^\r 



2k 



d9i 



,2 



Ayr 



k (jarpl"^ + \(3rs 



(11) 



k (\arp 
cot 9j 



A9x^ = 



2k [\arpf + \prs 
-f7[|rpp + |r^p -|-2|rp||r,|cos 
-i\l\ + l)(\af\rp\^ + m' 



2xks||7'p| sin( 



kzi 



(|arpp + |/3r. 



x(|/| + l)cot&,(|rp|'- 
2kzR (^arp\^ + \l5rs 



{12) 



(13) 



(14) 



where Vp^^ = |rp,s| exp[j0rp.J, ^ = 2Im[a*/3], x = 
2Re[Q!*,3], zji = kwl/2 is the Rayleigh length. When 
evaluating the energy term (E|E), we discard the cross- 
polarization terms (terms proportional to kx^/k or ky^/k) 
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FIG. 3. (Color online) The GH and IF shifts in TIR. The 
parameters are: incident angle 9i = 7r/4, refractive index of 
LHM riL = -1.515, refractive index of RHM ur = 1.515, 
vortex charge I — I, propagation distance z = zr, beam waist 
Wo = 20A, where A is the wavelength in vacuum. To make 
the shifts more noticeable, we amplify them by 20 times. In 
the first row, the incident beam is left-circularly polarized. 
The actual beam shifts of Fig. Ela) are {x^) — 1.129A and 
(yr) = -0.282A. The actual beam shifts of Fig. [S^b) are 
(x^) = -1.129A and (y^) = -0.282A. In the second row, 
the incident beam is mixed linearly polarized with x = 1- 
The actual beam shifts of Fig. [3jc) are (a;^ ) = 1.129A and 
{Vr) = 0.101 A. The actual beam shifts of Fig.g^d) are {x^} = 
-1.129A and {y?) = -O.IOIA. 




FIG. 4. (Color online) The spatial GH and IF shifts of trans- 
mitted vortex beams at air-LHM and air-RHM interfaces. 
The parameters are: incident angle 9i — tv/A, refractive index 
of LHM TiL ~ —1.515, refractive index of RHM ur — 1.515, 
polarization state a = 1/2, x = 1/2, propagation distance 
z = zr, beam waist wo ~ 20A. To make the shifts more no- 
ticeable, we amplify them by 200 times. In the first row, the 
vortex charge is i = 1. The actual beam shifts of Fig. Ufa) 
are Axf = 0.006A and Ay^ = 0.085A. Figure (Kb) has the 
same shifts with Fig. [3Ka). In the second row, the vortex 
charge is / — —1. The actual beam shifts of Fig. \M.c) are 
Axf = -0.006A and Ay^ = -0.033A. Figure (Kd) has the 
same shifts with Fig. [3](c). 




and Zr-component electric field E^^, since their average 
energy density is proportional to 9q. The first terms of 
Eqs. (fTT|) and are vortex-induced spatial GH and 
IF shifts. They were initially proposed by Bliokh [ISJ 
and Fedoseyev [i^, respectively. Experimental demon- 
strations were accomplished by Merano f46| and Das- 
gupta [11] at air-glass interface. Note that the vortex- 
induced spatial GH shift only exists in mixed linearly 
polarized beams (polarization vector is oblique to the in- 
cident plane), while the vortex-induced spatial IF shift 
occurs in arbitrary polarized state. The second term 
of Eq. Jlip coincides with the well-known Artmann for- 
mula [3]. The second term of Eq. (IT^ is the spin- 
dependent IF shift 0, \^ . In partial reflection and 
TIR, it turns into the Bliokh formula and Schillings 
formula respectively. Equations and are 
the angular shifts. In general, they are proportional to 
9q. For p-polarized beams, these angular shifts could be 
remarkably magnified through Brewster resonance [l7j 
(finally in the scale of 9o near the Brewster angle). In 
these cases, the cross-polarization terms are not negligi- 
ble when we calculate the energy intensity. 

A careful assessment of Eqs. pT|) - (fTil) indicates that. 



in partial reflection region, the spatial and angular GH 
and IF shifts at loss-free air-LHM and air-RHM inter- 
faces are identical. In TIR, shifts are different, how- 
ever. In this case, the vortex-induced shifts and angu- 
lar shifts are null. The phase 4>rp ^ has the same mag- 
nitude but the opposite sign for RHM and LHM, re- 
sulting in a negative GH shift [H [13 in LHM. The IF 
shifts for circularly polarized beams and mixed linearly 
polarized beams are —a cot6i[l + cos{4>r^ — (j)rp)]/k and 
— X cot 0i sin(0ra — 0r respectively. Here, we theo- 
retically predict that owing to the phase reversion, the 
IF shift in total reflection region would also be reversed 
at air-LHM interface when the incident beam is in mixed 
linearly polarized state. This prediction is illustrated in 
Fig. [21 With the refractive index nji = —ul = 1.515, 
we choose incident 9i = 7r/4 to avoid the deformation 
of reflected beams [11| . Figure [3l^a) and Fig. ^h) show 
that the IF shift in TIR remains unreversed when the in- 
cident beam is in circularly polarized state. Figure 
and Fig. ^d) suggest that the IF shift in TIR is reversed 
when the incident beam is in mixed linearly polarized 
state. 

Now, we begin to study the shifts of transmitted vortex 
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FIG. 5. (Color online) The spatial GH and IF shifts of 
transmitted vortex beams at air-LHM and air-RHM inter- 
faces. The parameters are: the refractive index of LHM 
riL = —1.515, the refractive index of RHM Tin = 1.515, po- 
larization state a = 1/2, x ~ V^- The spatial shifts are 
presented in the scale of A. 



FIG. 6. (Color online) The angular GH and IF shifts of 
transmitted vortex beams at air-LHM and air-RHM inter- 
faces. The parameters are: the refractive index of LHM 
riL ~ —1.515, the refractive index of RHM ur = 1.515, po- 
larization state a — 1/2, x = 1/2- The angular shifts are 
presented in the scale of Oq . 



in LHM. We substitute Eq. ^ into Eq. ^ and get the 
spatial and angular shifts of transmitted beams 



Xvl cot Oi (^\tp\'^ - \ts\' 



2k(^\atp\^ + \(3ts 

k[\atp\^ + \(3ts\^) 
k {jatpf + \l3tst 



(15) 



cot 9i 



'^XV\ts\\tp\ sin( 



2k [\atp\^ + \pts\ 
~a[\tp\^ + \ts\^ - 2Tj\tp\\ts\ cosi^t^ - <j>t^)] \,{16) 
Vm + l){\c^f\tp\%} + \l3f\ts\%^ 



kzR^ (|aip|' + |/3t,|' 



x{\i\ + i)coie,{\tp\' ~\t,i 

2kzR^ {\atp\^ + \l3ts\' 



(17) 



(18) 



where tp^s = |ip,s| exp[i(^tp J, zr^ = nkrfwl/2 and 
zr = nkwQ/2 is the Rayleigh length along xt and yt 
axis. These Rayleigh lengths of transmitted vortex beam 



in LHM are both negative owing to negative phase ve- 
locity Hi]. The first terms of Eqs. dm and ^ are 
vortex-induced spatial GH and IF shifts. Until now, no 
experiments have been reported on these vortex-induced 
GH and IF shifts. The vortex-induced GH shift of trans- 
mitted beam occurs when the incident beam is in mixed 
linearly polarized state. The vortex-induced IF shift, 
however, exists in any polarization state. The second 
term of Eq. (fTS)) is the phase-dependent GH shift. In 
lossy media, this part is not negligible. The second term 
of Eq. (dH) is the spin-dependent IF shift. If the LHM is 
loss- free, this term degenerates into Bliokh formula [11]. 
A special case is "total transmission", where n = —1. In 
this case, tp = tg = 1, we get {xt) = {yt} = 0. 

We closely examine Eqs. p5|) - ((T8)) and find that spa- 
tial GH and IF shifts of partial transmitted beams remain 
unreversed at air-LHM interface compared with air-RHM 
interface. The angular shifts, however, are reversed. Fig- 
ure m shows the spatial GH and IF shifts of transmitted 
beams at air-LHM and air-RHM interfaces. The beam 
profiles are drawn according to Eqs. (jA.4[) . (IA.7I) and 
(IA.8I) . We can easily find that the transmitted beam is 
stretched along xt axis, which is a pure geometrical phe- 
nomenon. We assume that the LHM is low-loss or loss- 
free. This assumption is reasonable considering state-of- 
the-art micro/nano manufacturing technology. For ex- 
ample, by incorporating gain media into fishnet struc- 
ture [49[ or exploiting second-order magnetic resonance 
of the fishnet structure 50], we can obtain low-loss, 3- 
D, polarization- independent LHM in the visible spectral 
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range. Based on this assumption, the phase-dependent 
term of spatial GH shift vanishes and the spatial GH shift 
only depends on vortex. This argument is confirmed by 
Fig. m^a) and Fig. Hl^c), where the vortex charge is ^ = 1 
and I = —1, respectively. They indicate that when the 
vortex charge is reversed, the spatial GH shift is also 
reversed. The spatial IF shift is both vortex and spin- 
dependent. Therefore, when the vortex charge is flipped, 
the absolute value of spatial IF shift is altered [but not 
flipped, see Fig.SJJa) and Fig.SJ^c)]. In Fig. [51 we demon- 
strates how incident angle impact the spatial GH and IF 
shifts. When the incident angle increases, both GH and 
IF shifts increase. The IF shift fails to converge when 6i 
approaches to 90°. We note that in this case, the energy 
transmission coefficient is zero. Hence, the spatial shifts 
of transmitted beams are meaningless. Figure [5] also sug- 
gests that both spatial GH and IF can be enhanced by 
raising the vortex charge. 

We proceed to analyze the angular GH and IF shifts. 
Figure [6] demonstrates how the incident angle affects the 
angular shifts. Figure [D^a) and Fig. IH^b) show the an- 
gular GH shift at air-LHM and air-RHM interfaces, re- 
spectively. They indicate that the angular GH shift is 
reversed in air-LHM. Figure [Gjc) and Fig. [6l^d) show 
the angular IF shift at air-LHM and air-RHM interfaces, 
respectively. They indicate that the angular IF shift is 
reversed in air-LHM. For angular IF shift in air-LHM, 
there exists a minimum point as incident increases from 
0° to 90°. These angular shifts increase proportionally 
as we raise the absolute value of vortex charge. 



IV. LINEAR AND ANGULAR MOMENTA 

In this section, we start to analyze the linear and angu- 
lar momenta of reflected and transmitted vortex beams. 
In general, the linear momentum density can be divided 
into orbital part and spin part [El'] , which reads 

p = Im[E* X (V)E] + ilm[V x (E* x E)], (19) 

where the first term is orbital momentum density p°, 
the second term is spin momentum density p**, E is the 
electric field in position space. In momentum space, we 
can write Eq. (IT9l) into 



p=\GE\^'k+{k,y-ky±){G*E;GEy-GE,G*E;). (20) 

Therefore, the momentum per unit length is 

P = (E|k|E) + {E,\k,\Ey)-y - {Ey\ky\E,)±, (21) 

where the first term is orbital momentum P*^, the sec- 
ond and third terms are spin momentum P'^. From the 
above equation, we find that Pz oc nk, which means Pz 
in LHM is opposite to Pz in RHM. For both paraxial and 
nonparaxial beams, the spin part makes no contribution 
to the linear momentum |44|. 15211: P'^ = 0. 



In position space, the angular momentum density de- 
notes j = r x p. In momentum space, we replace r with 
operator i9kj^ + zz and get the angular momentum 

J = (GE|(i9k^ +zz) X klGE) 

+ [{GEy\lidk^k, + dk^ky)\GE,) 

-{GKMdkJcc+dk^ky)\GEy}]z, (22) 

where the first term is orbital angular momentum J*^, the 
second and third terms are spin angular momentum J'^. 
Note that the transverse spin angular momenta are null, 
which means = Jy =0. We substitute the prop- 
agation operator G = exp \inkz — (/c^ + A;^)] into 
Eq. (I22p and get the orbital and spin angular momen- 
tum [ii, H il 

= (E|mfc9fe„ |E)x + (E| - inkdk^ |E)y 

+ {^^{ky^k^ - k^dk^)\E)i, (23) 

= [{Ey\i{dk^k^ + dk^ky)\E,) 

-{EMdk^K + dk,ky)\Ey)]i. (24) 

We find that the transverse orbital angular momenta 
and Jy are proportional to nfc, which means they have 
opposite directions in LHM compared with the counter- 
parts in RHM. Equations (|23p and ((24l) also indicate that 
the z-component orbital angular momentum and spin 
angular momentum have the same directions in LHM 
and RHM. Comparing Eq. HU]) with Eqs. ^ and 
we can easily get the following relations 



jO PO PO 



(25) 



where Aa; and Ay are spatial shifts, A9x and Ady are 
angular shifts. 

We first offer a qualitative explanation on why spatial 
shifts are unreversed in LHM but angular shifts are re- 
versed in LHM. Equation (|25p indicates that spatial shifts 
have no direct relations with the z-component orbital an- 
gular momentum and spin angular momentum . 
They are dependent on the transverse angular momenta 
J2 and z-component linear momentum P^ |^]. Since 
the directions of and P^ are both reversed in LHM, 
the spatial shifts Arj^ would remain unreversed as a re- 
sult. The angular shifts, however, solely depend on linear 
momentum. The reason why angular shifts are reversed 
in LHM is that the transverse linear momentum Pj is 
unreversed but the z-component linear momentum P^ is 
reversed. 

From now on, we will perform a quantitative analysis 
on the linear and angular momenta and confirm the mo- 
mentum conservation laws. Therefore, we assume that 
the LHM is loss-free and the amplitude reflection coef- 
ficients Tp^s and amplitude transmission coefficients ip s 
are real variables. By substituting Eqs. ([7]), and ([5]) 
into Eq. (j2ip . we obtain the linear momenta of incident. 
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reflected, and transmitted beams 



P -Zfcll^x +^fc^(r^-r^)cot0-y 
' 2zR 90, 2zR ^Jcotw.y, 

{\l\ + l)dT^ _ X(l^l + 1),,2 ,2^ ~ 

-[tp -t^) cot 9iyt 



2TfzR de, 

nkT 
^ zt, 



2T]Zf 



(26) 



where = Iq^'^pP 

|2 



l/^r^p is the energy reflection coeffi- 
cient, T = \atp\"' + \(3ts\'^ ■ The energy transmission coeffi- 
cient is Qt = \n\riT. Note that the linear momenta of in- 
cident, reflected, and transmitted beams are presented in 
three different coordinate systems Oi^r,t — Xi^r,tyi,r,tZi^r.t- 
We can verify that the linear momenta along x and y 
axes satisfy the conservation law [l^, [l^l 



P^. sin6'i — Pz^ smOi — P^__ cosf 



-\n\-q 



P — P 



■2(P,, cos^t-P,, sin^O, 



(27) 



For incident beam, there is no transverse linear momen- 
tum. To satisfy the linear momentum conservation law, 
transverse linear momenta Px^ , Py,. , Pxt i Pyt wo- 
duced. These transverse linear momenta are responsible 
for angular shifts. 

We proceed to analyze the angular momenta. By sub- 
stituting Eqs. ©, ©, and dH) into Eqs. ^ and (IM]) . 
we get the angular momenta of incident, reflected, and 
transmitted beams. 



3i ~ {I + a)i. 
J, = 



/ dQr cr , , o 

(r„ +rs)-^cot( 

2 ddi 2^ ^ ' 



'y(^p ~ ^s) cot6'iyr + (-^Qr + arprs)Zr, 



nl dT na , ^ ^ 



Xt 



l{l + Tf)T ^ atpts 



2ry2 



1] 



(28) 



where the spin angular momenta of incident, reflected, 
and transmitted beams are cr, arpTs, and atptsTj^^, re- 
spectively. For each individual photon, the angular mo- 
menta of incident, reflected, and transmitted beams are 
{I -h a)h, {-I + arprs/Qr)fi, + '7"^)/2 + crtpts/T]n, 
respectively. Though the orbital and spin momenta of 
reflected photon and transmitted photon depend on the 
absolute value refractive index, they are independent on 
the sign of refractive index. The z-component angular 
momenta satisfy the conservation law: 

Jz- cos 9i = —Jz^ cos 6i — Jx^ sinOi 

+ \n\7f{Jx,smet + Jz^cosOt). (29) 



Prom Eqs. (|28|) and (|291) . we can infer that there are 
two types of momentum conversions. The first type is 
spin-orbit conversion. The z-component spin angular 
momentum of incident beam a cos 9i converts into trans- 
verse angular momenta and J^, resulting in spin- 
dependent IF shifts. The second type is orbit-orbit con- 
version. The z-component orbital angular momentum 
of incident beam I cos 9i converts into transverse angu- 
lar momenta and , resulting in vortex-induced IF 
shifts. The vortex-induced GH shifts, though relate to 
Jy^ ^ , are not governed by angular momenta conservation 
law. 

In TIR, the linear and angular momenta of incident 
and reflected beams are 

P? — ^'Zj, Py- — fcz^, 

3i = {1 + a)zi, 

J,. = — cot 9i[a- + a C0s{(f>r^ — (t>rp ) + X Sni{(f)r^ — )] Xr 



2 



+ [ - ^ + cr cos(0r, - ) + X sin((/)r^ - cj)r^ )] , 

(30) 

We can easily verify that they fulfil the conservation law 



Pz, sin6'i = Pz^ sin 9i, 
Jz, cos Oi = —Jz^ cos t 



Jx^ sin ( 



(31) 



At this point, we would like to add three comments. 
First, it's worth noting that we adopted the Minkowski 
momentum in Eqs. ^7} and (|29p in this paper. Although 
the Abraham-Minkowski dilemma has been solved [53|- 
[55} . why Minkowski momentum is a proper form in this 
study, we believe, is an interesting problem worth further 
investigation. Second, two important papers [sgIIstI were 
published while our paper was being peer-reviewed. One 
paper |56|] adopted the real-space approach and revealed 
the role of longitudinal field Ez in IF shifts for the first 
time. We think additional analysis on the connection 
between longitudinal field and spin-orbit conversion can 
also be carried out in the momentum space. The other 
paper [s^l unambiguously separated the effects of beam 
shape and other parameters (such as polarization, the 
property of the interface) on GH and IF shifts, which has 
long escaped researchers' attention over the past years. 
But the discussions are confined to reflected vortex beam. 
Generalization to transmitted vortex beam, we think, re- 
mains quite challenging. Third, Owing to the close sim- 
ilarity between light beam and matter waves, scientists 
have found that electron beams can also possess orbital 
angular momentum by passing through a spiral phase 
plate [5^ or nanofabricated diffraction hologram [s^ . In 
this regard, we hopefully predict that vortex electron 
beam might also experience vortex-induced shifts in a 
potential well. By properly designing the potential well 
and taking advantage of quantum weak measurements 
or other measuring technology, we may even observe the 
vortex-induced beam shifts of electron beam in experi- 
ment. 
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CONCLUSIONS 



by inverse Fourier transformation 



In conclusion, we have derived the reflected and trans- 
mitted fields of vortex beam at air-LHM interface via an- 
gular spectrum method. By using this method, we have 
managed to get the formulas of spatial GH shifts, spa- 
tial IF shifts, angular GH shifts, and angular IF shifts. 
These formulas suggest that the spatial GH and IF shifts 
remain unreversed at air-LHM interface compared with 
air-RHM interface. By raising the vortex charge, the 
spatial shifts can be remarkably enhanced. In TIR, apart 
from reversed GH shift, we predict that the IF shift would 
also be reversed when the incident beam is in mixed 
linearly polarized state. The physical interpretation of 
these interesting phenomena lies in the reversed trans- 
verse angular momenta and reversed linear momenta. Al- 
though the spatial shifts have no direct relations with the 
z-component angular momentum (z-component angular 
momentum is not reversed in LHM), they are actually the 
outcomes of spin-orbit and orbit-orbit conversion. There- 
fore, the unreversed spatial shifts are indirect evidence of 
unreversed angular momentum of LHM. As for angular 
shifts, they are reversed at air-LHM interface and can be 
amplified by enhancing the vortex charge. This is direct 
evidence on the reversed linear momentum of LHM. Be- 
sides these qualitative analysis, we also offer concrete ex- 
pressions of the transverse linear and angular momenta, 
which explicitly reveal the physical picture of spin-orbit 
and orbit-orbit conversions. These momentum conver- 
sions are governed by z-component angular momentum 
conservation law. 



Appendix: ELECTRIC FIELDS IN POSITION 
SPACE 

In this appendix, we will give the full analytical ex- 
pressions of the electric fields of reflected and transmit- 
ted beams. The electric fields in position space are given 



u{x, y^^) = I dkxdkyu{kx, ky) exp[i(fc2;a; + kyy)] 



X exp 



inkz — ——iki + kl) 
2nk y 



(A.l) 



We first apply Eq. (jA.l[) to calculate the reflected flelds. 
Substituting Eq. (O into Eq. (|A.ip . we get the reflected 
beams 



9r„ d ^ , , ^ d 

P[rs + r.p) cot 0, — - 

dyr 



dOi dxr 
,9rc d 



d 



/^"^"^ +Oi{rs+ Tp) cot„, 

OOi OXr Oyr 



if d ^ d 

k \ OXr oyr 



Cikwo 



2(z_R + iZr) 



X exp 



kwQ Xr — i sgn[l]yr 

%/2 ZR+ iZr 



kjxj. + yj) 
2{zr + iZr) 



+ ikzr 



UrYr 

(A.2) 
(A.3) 



Note that the orbital angular momentum of reflected 
beam is —Ih per photon. Eqs. (jA.2[) and (jA.3|) are the 
full expressions of reflected electric flelds. 

We proceed to calculate the transmitted electric flelds. 
Substituting Eq. ([5]) into Eq. (lA.ll) . we get the transmit- 
ted fields 



E, 



atr 



f3ts 



q;77— -p— \- p [tp — rjts) cot 



d9i dxt 



dyt 
d 



dts d 

o9i dxt dyt 



if d n \ - 

—r atp- 1- pts-^ UtZt. 

nk \ dxt oyt. 



Ut^t 

utft 
(A.4) 



The formula of Ut is rather lengthy. If |^| is an even 
number, then the formula is 



Ut 



Ci\n\kwo 



1 



Wo 



V2. 



exp 



nk 
T 



E ^1 



—irj 


1 / — m 


'sgn[Z]" 









X iFi 



TO 1 



exp[-a62 - cA^]T 

1 -I- m 1 2 
— z — , cd 



* . -t- 




2/? 


. 




ZR 




IZtJ 


[1 + 1^1- 


TO 


r 


1 + TO 


2 


2 



inkzt 



—irj 




'sgn[lY 






[ \ 



X 1^1 



U| - TO 3 ,o 
2 ' 2' 



4:^/acbdexp[—ab^ — c(fi]T 



TO 3 ,n 

^ + Y'2'^' 





1^1 — TO 






TO" 






r 


n 






2 









(A.5) 
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if |/| is an odd number, then the formula is 
Ci\n\kwa 1 



Ut 



\/{zr.^ + izt){zR^ + izt) 

\i\-i 



Wo 



1 1'l 



V2. 



cxp 



nk 

T 



m=0,2,4- 





\l\—rn 











X iFi 



W 



1 + 



iFi 



2v^6exp[-a5^ - ccfjV 

1 + m 1 ,2 
-, ca 



zr^ + izt zr^ + izt 
l\ — m 



1 



+ inkzt 
1 + m" 



m— 1,3,5 



—irj 


\t\ — m 


'sgn[;]" 






[ Vc \ 



1 + 1^1 - 771 1 ,2 

—2 '2'"^ 



2 ' 2 

71 

2Vcdexp[-a62 - cd^jF 

771 3 ,n 

l + Y'2'^^ 



"14 


- 1^1 — 777 






m" 






r 


n 






2 









r 



where, C™ is the binomial coefficient, F is the Gamma 
function, iFi is the Kummer confluent hypergeometric 
function, a = 7;^7x;q/4 + izt/2nk, b = ixt/2a, c = Wq/A + 
izt/2nk, d = iyt/2c. Equations (jA.4[) - (|A.6p fully describe 
the transmitted vortex beams. Although Eqs. (jA.Sp and 
(|A.6p are cumbersome, we still manage to get the electric 
fields for several low order vortex beams 



izt{l - rf) 



+ 



nk{zR^ + izt){zR^ + izt) 

yxt iyt 

zr^ + izt zr^ + izt 
iztjl ~ 77^) 
nk{zR^ + izt){zR^ + izt) 



.1=1 



Wt 



cx 



zr^ + izt zr^ + izt 
iyt 



cx 



ZRy + IZt 

m 



ZR^ + IZt 

: 1 : — 

zr^ + izt zr^ + IZt 



(A.7) 
(A.8) 



(A.6) 



(A.9) 



(A.IO) 
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